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£NJ , Abstract 

A mixed Ising-Heisenberg spin system consisting of triangular XXZ-Heisenberg spin clusters assembled into a 
chain by alternating with Ising spins interacting to all three spins in the triangle is considered. The exact solution 
of the model is given in terms of the generalized decoration-iteration map and within the transfer-matrix technique. 

f*"*- ' Exact expressions for thermodynamic functions are derived. Ground state phase diagrams, thermodynamic and 

CNJ . magnetic properties of the system are examined. 

q ■ 1 Introduction 

Frustrated spin systems have been in the focus of intensive investigations during the last decades [T]. Belonging 
to a separate class of magnetic systems, they possess a number of fascinating features. Geometrically frustrated 
magnets are of special interest due to their unusual magnetic and thermodynamic properties. In this class of magnets 
the lattice topology precludes the spins from simultaneously minimizing all spin-spin-exchange interactions. The 
simplest example of frustrated spin cluster is given by the triangular arrangement of three spins with antiferromagnetic 
interaction between each pair. 

Numerous variants of lattice spin systems containing triangular plaquettes are, on one hand, widespread struc- 
tures of magnetic materials and provide, on the other hand, a theoretical prototype models for investigating the 
geometrically frustrated systems and their unusual features. Among the important classes of geometrically frustrated 
two-dimensional systems, it is worth mentioning the triangular antiferromagnets [2] and antiferromagnets on kagome 
lattice [3]. There are few one-dimensional lattices with triangular frustrated units which are famous for their ex- 
actly known dimerized ground states: the Majumdar-Ghosh model which is the special case of spin-1/2 chain with 
^y-s \ competing nearest-neighbor (NN) and next-nearest-neighbor (NNN) interaction [1] and the so-called sawtooth chain 

^ ' (A-chain), the system of corner-sharing triangles [5]. 

Although, properties of the ground state of the above-mentioned systems have been investigated analytically very 
well, only numerical calculations still remains reliable to gain an insight into the finite T thermodynamic properties. 
On the other hand, there is a large class of exactly solvable one dimensional "classical" lattice models which allow 
one to obtain an exact expressions for thermodynamic functions using different techniques [BJ. The simplest and at 
the same time key example is the transfer-matrix solution of the one-dimensional Ising model known since 1920s. 
Further developments of the transfer-matrix technique which is very powerful and rather simple for the arbitrary 
one dimensional lattice systems with discrete commuting variables turned out to be especially useful in the statistical 
mechanics of macromolecules|7] . However, in the theory of magnetism and strongly correlated systems one-dimensional 
systems with classical variables have rather "pedagogical" value than practical one. 

Nevertheless, it was established in a series of recent publications that considering the Ising counterparts of some 
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Heisenberg one-dimensional spin systems one can obtain, on one hand, the exact thermodynamic description of 
the problem and, on the other hand, the system which at least on the qualitative level exhibits the magnetic and 
thermodynamic behaviors very similar to those of the underlying quantum spin model [8]-|13j. For instance, the 
form of the magnetization curve does not acquire crucial changes when one replaces some or even all Heisenberg 
interactions in the chain with Ising ones. So, if the initial quantum spin chain exhibits a magnetization plateau at a 
certain value of the magnetization this feature will hold for the corresponding Ising or Ising Heisenberg chain, whereas 
quantitative characteristics of the plateau (terminal points, width, e.t.c.) can be different. This feature was reported 
in Ref. [8] for the spin-1/2 chain with bond alternating ferronragnetic-ferromagnetic-antiferromagnetic (F-F-AF) 
interaction, the model of 3CuCl2- 2dx (dx=l,4-dioxane) compound. In the quantum model describing the magnetic 
structure of 3CuCl2- 2dx, considered numerically in Ref. [U] the appearance of magnetization plateau at 1/3 of 
the saturation magnetization value was established. It is worthy to note that experiments have shown no plateau 
in 3CuCi2- 2dx which is caused by the insufficient small ratio of the antiferromagnetic and ferromagnetic coupling 
constants. The corresponding regime of magnetic behavior for F-F-AF Ising chain was also established in Ref. 8 This 



was the example of the intermediate magnetization plateaus in one-dimensional spin system which in recent decades 
received considerable attention both in theoretical and experimental aspects [I] [TB]. Hereafter, the Ising counterpart 
of this model was solved exactly in thermodynamic context in Ref. [5] demonstrating the same magnetization plateau 
at m — 1/3. Almost the same program but with different techniques was performed in Ref. [9] for another bond 
alternating chain, the spin-1/2 F-F-AF-AF chain which is regarded as the model of magnetic structure of Cu(3- 
Clpy) 2 (N 3 )2, where 3-Clpy indicates 3-chloropyridine. The authors considered a simplified Ising-Heisenberg variant of 
the F-F-AF-AF chain within the method of exact mapping transformation (decoration-iteration transformation) [T7] . 
Comparing their results with the experimental and numerical data they found not only qualitative but also quantitative 
agreement as well. It is also worthy to mention the works on exact solutions of mixed spin-(l/2, 1,3/2) and spin-1/2, 
the so called diamond chains with the decoration-iteration transformation technique |10j . which revealed series of 
magnetization plateaus in magnetization process of the system. Also, the magnetization plateaus in one-dimensional 
spin-1, spin-3/2 and spin-2 Ising chains with single-ion anisotropy have been investigated in Ref. pj] showing a good 
agreement with experimental data obtained for spin-1 Ni-compounds. Magnetization plateaus in the Ising limit of the 
multiple-spin-exchange model on the four-spin cyclic interaction have been considered in Ref. [12] . 

All these results, unexpected at the first glance, offer the challenge to develop a formalism which, on one hand, 
leads to the exact analytical treatment for thermodynamic functions of the one-dimensional spin systems and, on the 
other hand, could be relevant in describing and understanding experimental data for real materials. As mentioned 
above, the idea is very simple, just to replace some (or even all) Heisenberg-type quantum interaction with the 
Ising ones. Of course, there are many very important features of quantum exchange interaction which are crucial in 
understanding many phenomena in magnetic systems and which cannot be neglected. But the results of Ref. [8]- 
[12] indicate the existence of at least a few quantum spin chains (actually the chains with alternating ferromagnetic 
and antiferromagnetic couplings), the thermodynamic properties of which do not undergo sizable changes under the 
replacement of some interaction with more simple Ising ones, allowing an exact solution. In addition to that, the spin 
systems considered here may serve as more or less realistic model for the arrays of single-molecule magnets (SMMs) 
[18] with inter-molecular couplings, engineered by the methods of the so-called supra-molecular chemistry [19]. SMMs 
recently receive much attention as they are very convenient objects for studying the phenomena which are critical at 
the nanoscale, such as quantum tunneling of magnetization, decoherence, and Berry phases. One of the well-studied 
specific examples of SMMs, the equilateral spin triangle CU3 [20] . assembled to the one-dimensional arrays with very 
weak inter-molecular coupling could be the system which approximately can be described by the model, considered 
here. 

In this paper we consider a further development of the methods applied in the Ref. [8]-[13j. for the chain with 
frustrated triangular quantum XXZ-Heisenberg plaquettes. This chain was introduced in Ref. [14] where the ground 
state properties and some numerical results for purely Heisenberg case were obtained. Within the transfer-matrix 
technique we obtain exact thermodynamic solution of the chain consisting of such spin triangles alternating with the 
single Ising spins which are coupled by the interaction of Ising type with all six spins of the sites of two neighbor 
triangular plaquettes. The paper is organized as follows. In the Second section we formulate the model and present 
its solution by the decoration-iteration transformation technique. In the Third section the transfer-matrix solution of 
the system is given. In the Next section, using the transfer-matrix formalism, the exact thermodynamics of the system 
is discussed. In particular the T = phase diagram and plots of magnetization vs. external magnetic field, specific 
heat and magnetic susceptibility are presented. The Appendix contains several technical points of the calculations. 

2 The model and its solution with the aid of the generalized decoration- 
iteration transformation. 

We begin with the hybrid Ising-Heisenberg spin system consisting of the triangular clusters of Heisenberg S — 1/2 
spins interacting to each other with coupling constant J and axial anisotropy A. These triangles are assembled into a 
chain by alternating with single sites. Each spin at the single site is coupled to all six spins situated at the corners of 
its two adjacent triangles (see Figfl]). The interactions of spins at single sites and spins belonging to triangles are all 
supposed to be of Ising type (the interaction includes only z-components of the spins). The corresponding Hamiltonian 
is suitable to be presented as a sum over the plaquette Hamiltonians each one containing one triangle and its two 



surrounding single sites: 
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Here by a we denote the ^-components of the spins at the single sites which, thus, can be identified with the Ising 
variables taking values ±1, K is Ising coupling between triangle spins and single site spins, H\ and H 2 stand for the 
couplings of Ising and Heisenberg spins to the external magnetic field pointing in z direction. For convenience we 
normalize the Heisenberg spin operators in such a way which provides the eigenvalues of S z to take ±1 values, which 
means that our S a are just the Pauli matrices without spin magnitude multiplier. (In order to recover conventional 
Heisenberg Hamiltonian parameters one should replace in all our formulas J with J/4 and H with H/2) Thus, they 
obey the commutation relations: 

[Sta^r b ] = 8 ij d ab 4Sf a , (2) 

and act on the basic states in the following way 

■S*|T> = |T>, 5«|i) = -|i), (3) 

s+|T) = o, s+||) = 2|T) 

5-|T> = 2||), 5-||) = 0. 
In order to describe thermodynamics of the system one need to calculate the partition function 

Z = £Tr (s) e-^-"*, (4) 

where the sum is going over all possible configurations of the Ising spins and Tr(g) denotes the trace over all Heisenberg 
operators S and j3 as usually is the inverse temperature. One can easily see that the Hamiltonians corresponding to 
different plaquettes do commute which allows one to expand the exponent in the partition function obtaining partial 
factorization of the partition function 

N 
(a) t=l 

here Tr^ stands for the trace over the state of ith triangle. Now we can utilize the so-called decoration-iteration 
transformation [10J which allows us to represent the trace over all states of single triangle interacting with two Ising 
spin as the elements of the transfer-matrix for ordinary one-dimensional Ising model with a certain renormalization of 
its parameters: 

8 

here A n (<7j, <Ji+\) are eight eigenvalues of the TLi which depend of the values of its neighbor Ising spins. The parameters 
of emergent "one-dimensional Ising" model are expressed by the parameters of the initial model in the following way: 

A=(Z + Z_Zt) 1, \ /?fl=ilog(^), /3 # = /3Fl + ~log(J±), (7) 

where 
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are the partition functions for single triangle calculated for certain values of the its neighbor Ising spins. The explicit 
form of these functions can be found in the Appendix. Thus, according to Eqs.(TJ}-([5]) one can obtain an exact solution 
for the system under consideration in the form corresponding to the well known solution of one dimensional Ising 
model as follows: 
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Z(J3; J, K, A, H U H 2 ) = £) [] J 4 e /M* i a 4+1 +/9f(«r 1 +«r< +1 ) = A N Z (fc R, H), (9) 
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The partition function of the one-dimensional Ising model with periodic boundary conditions in thermodynamic limit 
is® 

Z (/3; R, H) = e pNR (cosh((3H) + ^sinh 2 (0H) + e~^A . (10) 

The relation (|9|) establishes a connection between thermodynamic quantities of the system under consideration and 
those of one dimensional Ising model. So, in the thermodynamic limit for free energy per spin we get 

(We remind that by N we denote the number of the lattice blocks each one containing five spins, so under the periodic 
boundary conditions the total number of spins is 4N.) Here by /o we denote the corresponding free energy per spin 
for ordinary Ising chain. Now we can utilize the common thermodynamic relations to obtain expressions for various 
thermodynamic quantities we are going to investigate. As usual for specific heat, entropy, total magnetization and 
susceptibility per spin one has 

<>-<&).< *--{%).- "--(!&■ *=(S) T 

The thermal averages of various microscopic variables can also be expressed in term of corresponding quantities for 
the Ising model. The sublattice magnetization and various nearest neighbor correlation functions can be expressed by 
the sublattice magnetization and nearest neighbor correlation functions of ordinary Ising chain given by 
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For sublattice and total magnetizations for the model under consideration one obtains: 
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where the corresponding functions W a , a — ±, can be found in Appendix. In the same manner for nearest neighbor 
correlation functions we obtain 

Caa = (<JiCFi+l) = C , (16) 

Cs, = (S ia ai) = — - [ M ) TH = -12 ( y + + 2V - M ° + V + C ^ > 

CIs = <Sf S£) = | (|j) " 2A 47 5 = ^ (^+ - F + + 2 (^ u ° F °) M ° + i^ u - F -) c °) . « ^ & - 
3 Solution within the transfer-matrix technique 

We will use another technique for describing the thermodynamics of the system under consideration. Namely, having 
the analytical expression for the single quantum triangle (see the Appendix) , one can represent the partition function 
of the chain under consideration given by Eq. (|5|) , in the form which mimics the partition function of the classical 
chain with two state variables at each site (the "Ising chain" with arbitrary Boltzmann weights): 

N 

Z = E II e^^+'^Z [fn, a i+1 ) = TrT w , (17) 
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where Z (tr^, o"i+i) is calculated in the Appendix and T is the following transfer-matrix : 



e^Z + Z 

Z Q e-^Z. 



(18) 



The eigenvalues of T are 



A± = ~ leP Hl Z+ + e-^Z- ± \J(e^Z + ~ e-^Z^) 2 + 4Z 2 J . (19) 

Thus, the free energy per one spin in the thermodynamic limit when only maximal eigenvalue survives is 
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For sublattice and total magnetization one obtains 
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where the functions X + , A_ , Xq are defined in the Appendix. If one assume H i = Hi = H which means the equality 
of the ^-factors for S- and er-spins the expression for the total magnetization per spin takes the following form: 
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where H2 — H should be put in all X a and Z a functions. 



4 Thermodynamics and magnetic properties. 

In this section we examine the thermodynamic and magnetic properties of the model under consideration. Having 
the exact explicit expressions for the magnetization, specific heat, entropy and nearest neighbor correlation functions 
one can easily obtain their plot vs. temperature or external magnetic field, despite of their complicated form. Also, 
analyzing the low-temperature region one can obtain the ground state properties and corresponding phase diagrams. 
In that follows, we restrict ourselves with the case of all antiferromagnetic couplings, J > 0, K > and A > 0. 

4.1 Zero temperature ground state phase diagrams. 

Let us analyze the possible ground state of the system at T — and for positive values of all model parameters J, K, A. 
According to the possible states of three s = 1/2 spins one can distinguish the following spin configurations on the 
chain under consideration. First of all, the ground states in the absence of magnetic field depending on the values of 
model parameters can be either ferrimagnetic (F) or antiferromagnetic (AF), 

TV 

|^>=ni3/2,3/2>,(g)||> l , (23) 

i=l 
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|AF)=n 11/2,1/2},® | 4),, 



where |Z, m)i stands for the spin state of triangle at i-th plaquette with total spin equal to I and S z = m and | T). | T) 
are the standard "up" and "down" states of the Ising spins. Expanding the states of quantum triangles by the single 
spin basis one finds 

|3/2,3/2> = | TTT), (24) 

|l/2, l/2)fl,i = -j= (| TU) + e ±2 ^| TIT) + <^| TTT)) , 

where the choice of the sings of the coefficients in 1 1/2, 1/2) state is connected with the its chirality. For the inter- 
actions considered here these two states are degenerated by energy and where is no need to distinguish them in our 
consideration, so we will omit R 1 L indices. Thus, here the manifestation of frustration consists in this degeneracy. 
Every triangle can spontaneously pass from 1 1/2, 1/2) ^ to 1 1/2, 1/2)^ and vice versa. Therefore, due to the frustration 
the system possesses non zero entropy at T — 0. It is worth mentioning, that considering the Ising limit of the system, 
A = 0, one can obtain essential enhancement of frustration at J — K. Indeed, if one takes only one plaquette which 
consists of five Ising spins with nine uniform couplings between them, the ground state will be disordered one with 
three spins pointing up and the remaining two pointing down or vise versa. So, the number of possible degenerated 
configuration for one plaquette and for fixed orientation of the total magnetization will be 10. The energies per one 
plaquette and magnetizations corresponding to \F) and | AF) configurations are: 

e F = 3r) - 2h - 6, M F = 1/2 (25) 

e AF = -v(l + 2A) - 2, M AF = 

where the following dimensionless parameters are introduced, e = E/NK, i] = J/K, h = H/K. Thus, the ground state 
of the system at T — and h — is ferrimagnetic when r\ < 2+^ and antiferromagnetic otherwise. The corresponding 
phase diagram is presented in Fig. ([2]). The appearance of the external magnetic field makes another two ground 
states possible. Namely, the states with total magnetization M = 1/2 resulting from AF-state by the flip of all Ising 
spins, which actually is a ferrimagnetic state differing from the (F) and completely polarized saturated phase: 

N 



|F) = ni 1 /2,V2) i 0|T) i , (26) 
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|5> = ni3/2,3/2) 4 ®|T)i, 



with 

e F = -77(1 + 2A) -2/i + 2, Mp = 1/2, (27) 

£<? = Srj - Ah + 6, Ms = l. 
The corresponding T = ground state phase diagram in the (h, rj) plane for the case A = 1 is presented in Fig. ([3]). 



4.2 Magnetization processes and susceptibility 

According to the phase diagram displayed in Fig. $5§ one can expect three different kinds of magnetic behavior 
depending on the relation between parameters 77 and A. The following sequences of transitions occur when the 
magnitude of external magnetic field is increasing. For the values of 77 and A belonging to the first region given 
by 77 < 2+A' ^ ne system undergoes only one transitions from \F) to \S) taking place at h = 6. The next range of 
parameters values is characterized by ^xa < V — 2+A' Here one can occur two consecutive transitions from \AF) to 
\F) at ft, = ?7(2 + A) — 2 and from \F) to \S) at ft = 6. And, finally, when 77 is greater than -^K ^ ne system undergoes the 
following transitions \AF) — ► \F) and \F) — > IS*) at ft. = 2 and ft = 77(2 + A) + 2 respectively. The corresponding plots 
of the magnetization processes for finite temperatures are shown in Fig.Q-©. At zero temperature all transitions 
are obviously jump-like and the corresponding magnetization curves are perfectly step-like. However, at arbitrary 
finite temperatures all transitions are smeared out within some interval of magnetic field. In Fig. (0J the plots of M 
vs H/K are presented for the r\ = 0.5 and A = 1 for several temperatures. One can see the magnetization plateau at 
M = 1/2 which corresponds to the zero temperature stability of \F) state in the interval of ft from to 6. The plot for 
T ' I K — 0.01 is very close to this picture, whereas one can see essential shrinkage of the plateau width with increasing 
temperature. Surely, for large enough temperatures the magnetization curve acquires the form corresponding to the 
free spins (Langevin curve) with monotonic increase between and 1 values. Generally speaking, at the small values 
of the ration 77 the system under consideration is similar to the mixed antifcrromagnetic Ising chain with alternating 
s — 1/2 and s — 3/2 spins, because the strong interaction between the triangle spins makes their behavior very similar 
to that of one spin-3/2. So, the region 77 < -^TE can ^ e considered as the spin-3/2-spin-l/2 mixed chain regime. The 
typical plots for the second range of the 77 is presented in Fig. (fSJ). Here two magnetization jumps corresponding to the 
transitions \AF) — > \F) and \F) — > IS") occur at T — in the magnetization curves for the system under consideration 
for rj = 1 and A = 1. The corresponding plateau at M — 1/2 as in the previous case appears due to stability of the 
\F) state between \AF) and \S). The width of the plateau is 8 — 7/(2 + A). The plots for different finite temperatures 
from Fig. [5] exhibit the gradual shrinking down of the plateau with increasing temperature. The same plateau at 
M = 1/2 but with another physical content one can see in the magnetization curves for the third region of values of 77. 
In the Fig. ([6]) one can see the plots of M vs. H/K for 7; = 1.5 and A = 1. All properties mentioned above concerning 
the temperature effect holds in this case as well. However, the plateau at M — 1/2 now corresponds to the \F) state 
of the chain. The two terminal points of the plateau correspond to the zero temperature transition from \AF) to \F) 
and from \F) to \S) respectively. The position of the left end of the plateau (ft = 2) is unaffected by the any changes 
of interaction parameters provided the later stay within the 77 > ^.( The same phenomenon for the right end of 
plateau at ft = 6 take place in the two other regimes of behavior of the system under consideration.) The width of 
this plateau at T = equal to 77(2 + A). 

The plots of thermal dependence of magnetic susceptibility times temperature T\ per one spin are displayed in 
the Figs. ([TJ)- (J9j - First of all, zero-field susceptibility diverges at T = which is a consequence of the fact that 
T = is the critical point at which the system possess an ideal long range spin order, which is destroyed by any finite 
temperature. In Fig. ([7]) the plots of T\ vs. T for the parameters first region are presented for several values of 
external magnetic field. Here no divergence is observed at T = because of energy gap opened by magnetic field. A 
rather sharp peak appears in the curve at extremely weak field value. This peak corresponds to the thermal instability 
at the transition from uncorrelated disordered state, which is the ground state at any finite temperature and H = 0, to 
the ferrimagnetic state \F). At the end of the peak the tiny round minimum is observed. With increasing of the field 
strength the low temperature peak is smoothed out and the curve become monotonically increasing, which is general 
feature of antiferromagnetically coupled spin systems. At H/K = 5.9, i.e. in the vicinity of the transition from \F) 
to \S) the Tx vs. T curve exhibits a very rapid increase at low temperatures with further narrow plateau formation 
which is then followed by the slowly and monotonically increasing region. The inset of Fig. ([7]) shows the small round 
minimum which accompanies the transition from narrow plateau to increasing monotonic behavior of susceptibility. 
Very similar picture of susceptibility thermal dependence was obtained for the F-F-AF-AF Ising-Heisenberg chain in 
Rcf. 9 . In Fig. (8) the plots of T\ vs. T for 77 — 1, A — 1 are presented. Here, a little smaller than in the previous 
case, low temperature peaks appear at the values of H/K corresponding to the transitions from \AF) to \F) and from 
\F) to \S). Similar picture one can observe for the Tx thermal behavior for the system under consideration when the 
parameters belong to the third region, 77 > ^4-^- The corresponding plots are presented in Fig. ©. 

4.3 Specific heat 

As mentioned above, the technique developed in the paper allows one to obtain analytical exact expressions for all 
thermodynamic functions of the system. Thought, the forms of some of them are extremely cumbersome, one can easily 



obtain their plots. In Figs. (|T0|) -(|T2 |) the plots of thermal behavior of the specific heat of the system under consideration 
are presented. Generally, specific heat exhibits two-peak structure, with one sharp low temperature, which is almost 
unaffected by the value of the magnetic field, and another one with broaden peak, which undergoes considerable 
changes under increasing magnetic field. Fig. (|I0[) demonstrates the corresponding plots for r\ = 0.5, A = 1. The 
second peak is well pronounced here and is strongly field dependant. Thus, one can associate it with Schottky-type 
anomaly. With the increase in the external magnetic field strength the peak moves to higher temperature region and, 
at the same time, drops in magnitude. At the values of field strength close to the transition to saturated state at 
T = the second broad peak begins to enhance again for a while, however with the further increase in the field it 
gradually gets smoother. In general, one can observe almost the same features in the thermal behavior of the specific 
heat times temperature per spin for the region ^ta — V < 2TA (FiglHI- The essential difference with the previous 
case is completely merging of sharp low-temperature peak and the Schottky-type broad peak at the values of field 
strength corresponding to the plateau in the magnetization curve at T = 0. Fig. (fT2j) shows typical plots of T\ vs. T 
for r\ = 1.5, i. e. for the third region of model parameters, 77 > tta- Here at low field strength immediate after the 
first sharp peak a narrow almost horizontal region of the curve appears. Increasing the field strength one can observe 
the enhancement of the slope of the corresponding region accompanying the second broad peak changes described 
above. 

5 Concluding remarks 

In this paper we considered exactly solvable model of one dimensional spin system with Heisenberg XX Z triangular 
spin clusters alternating with single Ising spins. Two different approaches have been discussed, the decoration- 
iteration transformation [9] QUI [UJ and the transfer-matrix direct formalism. Both approaches lead to the possibility 
of obtaining exact analytical expressions for the thermodynamic functions of the system. The system considered 
here, as well as the other systems considered earlier in Ref. [5]-[T2] at first glance is only of academic interest. 
However, the strict indications of the relevance of the approaches based on the simplification of the spin exchange 
interaction scheme in one-dimensional magnetic materials have been given in Ref. [H [9] [11] . Namely, considering 
the spin exchange Hamiltonians for various materials with one-dimensional exchange structure in most cases one can 
gain insight in description of their properties and understanding of their thermodynamic behavior only within the 
complicated and laborious numerical calculations demanding some time extremely powerful computing facilities, such 
as supercomputers. On the other hand, at least for some class of systems, the simplification of the model based 
on the replacement of all or just some of the exchange Heisenberg interactions with the interactions of Ising type, 
does not affect the prominent properties of their thermodynamics. This fact makes the simplified exactly solvable 
(in the "thermodynamic" sense) counterpart of the one-dimensional models which describes the magnetism of real 
materials very promising. Especially this approach can be successful in the series of bond alternating chain with 
spatially repeated sequence of ferromagnetic and antiferromagnetic interactions as it has been established for F-F-AF 
[8] and F-F-AF-AF[9] chains. Apparently, replacing the ferromagnetic bounds by Ising ones is especially "safe" in 
that sense, because the former ones favor parallel orientation of the spins, which is given by a separable vector of 
state. Thus, one can lose only a small amount of information of the properties of that bound by regarding it as the 
Ising one. Antiferromagnetically coupled spins are not so " harmless" , because they can be in the superposition and 
entangled states which cannot be properly translated into the Ising language. Another essential difference between 
Heisenberg spin chains and their simplifying counterparts is the spin-wave properties. Generally speaking, spin wave 
cannot propagate through Ising spins, so only localized excitations are relevant in that case. Further applications of 
the methods discussed in the present paper in other one-dimensional and quasi-one-dimensional Heisenberg models, 
especially in the models of novel magnetic materials, as well as a comparison of the results with the experimental data, 
can play a very important role in the understanding of complicated magnetic materials and their thermodynamic 
properties and can draw the frameworks of applicability of the approach offered in this paper. 
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7 Appendix 

Here we derive the exact expressions for the functions Z+, Z-, Zq from Eq. IJ5J). For these purposes one should obtain 
the eigenvalues of the matrix, corresponding to the Hamiltonian of single triangle. So, we define a partition function 
for one single triangle 

Z A (/3;J,A,H)^Tr e -^« A ^ = ^(^|e-^« A ^|^), (28) 

where \ip n ) is an arbitrary orthogonal normalized system of states. Then, one can easily see that the Hamiltonian 
corresponding to i-th plaquette from Eq. JT]) has the same form as the Hamiltonian of single isolated triangle if one 
introduces instead of magnetic field H an "effective field" depending on the values of o~i and <7j+i: 

H = H 2 -K{a t + n l+l ). (29) 

The term corresponding to the interaction of o~i and o^+i with magnetic field Hi adds — Q(ai+o~i+i) to each eigenvalue 
of H a ■ It is easy to see that Ha is diagonal in the symmetry-adapted basis the higher weight state of which is given by 
Eq. (|25|) . The rest 5 states can be easily found by successive action of the lowering operator S^ ot — ^(S^ + S 2 + S$) 
on |3/2, 3/2) and |l/2,l/2) i]fl : 

(30) 



|3/2,l/2) = 


^ (1 TU> + 1 UT> + 


litt)), 


13/2,- 


-1/2) 


= ^(iui) + im) 


+ IUT», 


13/2,- 


-3/2) 


= mi), 




11/2,- 


-l/2) fi , L = -i=(mT)+e ± 


^|TU>- 



Taking into account Eq. (|29[) the eight eigenvalues of the Hi are 

Ai )2 (o-i,o- i+ i) = 3(J± H) T 3K{a t + a l+1 ), (31) 

\z,i{o- l ,a l+ i) = A 5:6 (o-i,cr i+ i) = -(1 + 2A)J±H=fK(ai +a i+ i), 
A 7) 8(o-i, o-i+i) = -(1 - 4A)J ±Ft K{<ii + C7i + i). 

Thus, one have 
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Z(o-i, o-i+0 = ^ e^"^^ = Si cosh(/3( J ff 2 - ^(a, + CTi+1 ))) + B 2 cosh(/33(ff 2 - tffa + <7 l+1 ))), (32) 

i=l 

where 

Bl = 2 ( e /3(l-4A)./ + 2e /3(l+2A)J^ ( ^ = 2e -/33J^ (33) 



Now, we can define the functions appearing in Eqs. (jl5[) - (Tl6j) : 

Z+ = B x cosh(P(H 2 - 2K)) + B 2 cosh(/33(iJ 2 - 2K)), (34) 

Z_ = Si cosh(f3(H 2 + 2K)) + B 2 cosh(/33(ff 2 + 2K)), 
Z = Bi cosh((3H 2 ) + B 2 cosh(/33iJ 2 ), 

X »= T (^r) =B 1 smh(j3(H 2 -a2K)) + 3B 2 smh{[33(H 2 -a2K)) 1 a = ±,0, 



dH 2/T 
(dZ±\ 

T 



Y±=T(-^j -- t2 (Si smh(P(H 2 T 2#)) + 3B 2 sinh(/33(ff 2 T 210)) , 



<2« = T i^TT ) = B 3 cosh((3(H 2 - a2K)), a = ±, 0, 



V<9A 



r 



where 



Pa=T \~§T) = B 4 cos HP( H 2~a2K)) + B 5 cosh{P3{H 2 ~a2K)), a = ±,0, 

f± ~p; + z: ±2 t > Fa ~z;-z:> 



B 3 = 16Je /3 ( 1 - A ) J sinh(/33AJ) (35) 

B 4 = 2 ((1 - 4A)e^ 1 - 4A > J + 2(1 + 2A)e^ 1+2A ) J ) , 
B 5 = -6e^ 3J . 
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Figure 1: The schematic picture of the Ising-Heisenberg chain with triangular plaquettes alternating with single sites. 
The fillcd(cmpty) circles indicate Heisenbcrg(Ising) spins. 
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Figure 2: T = 0, H = ground-state phase diagram in the A — r\ plane. The boundary between ferrimagnetic(F) and 
antiferromagnetic(AF) ground states are given by r\ = 2+S 
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Figure 3: T — ground-state phase diagram in the rj — h plane for A = 1. The boundary between ferrimagnetic(F) 
and antiferromagnetic(AF) phases is h = f](2 + A) — 2; between second ferrimagnetic phase(F) and saturates state(S) 
\sh = 77(2 + A) + 2. 
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Figure 4: Total magnetization M as a function of dimensionless external magnetic field h for r\ — 0.5, A = 1 and three 
different dimensionless temperatures T/K = 0.01,0.5 and 1.2. 
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Figure 5: Total magnetization M as a function of dimensionless external magnetic field h for rj = 1, A = 1 and three 
different dimensionless temperatures T/K = 0.01,0.5 and 1.2.. 
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Figure 6: Total magnetization M as a function of dimensionless external magnetic field h for r\ = 1.5, A = 1 and three 
different dimensionless temperatures T/if = 0.01,0.5 and 1.2. 
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Figure 7: The temperature dependence of the T% product for several values of field strength and rj = 0.5, A = 1 Inset 
displays the round minimum for H/K — 3 on an enlarge scale . 
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Figure 8: The temperature dependence of the T\ product for several values of field strength and ry = 1, A = 1 
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Figure 9: The temperature dependence of the T\ product for several values of field strength and r\ = 1.5, A = 1 



O 



H/K=0.1 
H/K=3 
H/K=6 
H/K=8 




Figure 10: Specific head per one spin as a function of the dimensionless temperature T/K given in K units for 
ry = 0.5, A = 1 for various values of h. 
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Figure 11: Specific head per one spin as a function of the dimcnsionless temperature T/K given in K units for 
r\ = 1, A = 1 for various values of h. 
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Figure 12: Specific head per one spin as a function of the dimcnsionless temperature T/K given in K units for 
rj = 1.5, A = 1 for various values of h. 
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